Sound Velocity Anomaly at the Mott Transition: 
application to organic conductors and V2O3 
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Close to the Mott transition, lattice degrees of freedom react to the softening of electron degrees 
of freedom. This results in a change of lattice spacing, a diverging compressibility and a critical 
anomaly of the sound velocity. These effects are investigated within a simple model, in the framework 
of dynamical mean-field theory. The results compare favorably to recent experiments on the layered 
organic conductor k-(BEDT-TTF)2Cu[N(CN)2]C1 . We predict that effects of a similar magnitude 
are expected for V2O3, despite the much larger value of the elastic modulus of this material. 



X 



The Mott transition, which is the metal-insulator 
transition (MIT) induced by electron-electron interac- 
tions, has been investigated theoretically and experi- 
mentally for many years 0, 0. Some materials are 
poised very close to the Mott transition, which can 
therefore be induced by varying pressure, tempera- 
ture, or chemical composition. This is the case of 
(Vi_ a; Cr a ;)203 , and of the family of layered molec- 
ular crystals k-(BEDT-TTF) 2 X where X is an anion 
(e.g.,X=I 3 ,Cu[N(CN) 2 ]Cl,Cu(SCN) 2 ). In these com- 
pounds, one observes a pressure induced, finite tempera- 
ture, first order phase transition. Pressure increases the 
bandwidth, reducing the relative interaction strength. 
The first-order transition line ends at a second-order crit- 
ical endpoint (P c , T c ). The critical behaviour at this end- 
point has been recently the subject of remarkable exper- 
imental investigations 00. 

On the theory side, our understanding of the Mott 
transition has benefitted from the development of dy- 
namical mean-field theory (DMFT) 00. In this the- 
ory, electronic degrees of freedom are the driving force 
of the transition, and the critical endpoint is associated 
with a diverging electronic response function Xei (defined 
below) ■ An analogy exists with the liquid-gas tran- 
sition. The insulating phase is a low-density gas of neu- 
tral bound pairs of doubly occupied and empty sites; the 
metal is a high-density liquid of unbound doubly- occu- 
pied and empty sites which therefore conduct. The scalar 
order parameter is associated with the low-energy spec- 
tral weight. 

In real materials however, lattice degrees of freedom do 
play a role at the Mott transition. This is expected on 
a physical basis: in the metallic phase the itinerant elec- 
trons participate more in the cohesion of the solid than 
in the insulating phase where they are localized. As a 
result, the lattice spacing increases when going from the 
metal to the insulator. A discontinuous change of the lat- 
tice spacing through the first-order metal-insulator tran- 
sition line is indeed observed in (V]_ x Cr x )203 0. Ra- 



man scattering experiments |l(J find that in the metallic 
state the frequency of certain phonons associated with 
the BEDT-TTF molecules has a non-monotonic tem- 
perature dependence below 200 K. The effect of the 
Mott transition on optical (Einstein) phonons was stud- 
ied theoretically in |llj. Acoustic experiments find 
an anomaly in the sound velocity of the organic mate- 
rials as a function of temperature, with a particularly 
dramatic reduction recently observed 0] for k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Clnear the Mott critical endpoint at 
T c ~ 40 K. 

In this paper, we propose a simple theory of the effects 
connected with lattice expansion through the Mott tran- 
sition, and the associated divergence of the compressibil- 
ity. We address in particular the critical anomaly of the 
sound velocity observed in acoustic experiments. For this 
purpose, both the electronic degrees of freedom and the 
ionic positions must be retained in a model description. 
We adopt the simplest possible framework, previously 
introduced in Refja, namely the compressible Hubbard 
model (see also [lj]), with the electron correlations be- 
ing treated within DMFT. Our results compare favor- 
ably to the recent acoustic experiments on the layered or- 
ganic conductor k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 [referred 
to in the rest of this paper simply as the organic- 
conductor(OC).] Furthermore, we show that effects of a 
similar magnitude are expected for 2o3, despite the much 
larger value of the elastic modulus of this material. 

In the following, we assume that the dependence of the 
free energy on the unit-cell volume v (or, rather on the 
induced strain, see below) can be written as: 

F = F Q - P (v -v ) + l B {V - Vq)2 + F el [D{v)] (1) 

The last term F e i, is the contribution of the electronic 
degrees of freedom which are active through the transi- 
tion (e.g the d-shell for 2o3). Specifically we take F e i to 
be the free energy of a single band, half-filled Hubbard 
model with a half-bandwidth D{v) depending on the 
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FIG. 1: The variation of the electronic response function Xei 
with the bandwidth D for various temperatures. The param- 
eters are U/D = 2.492 and T/D = 0.060, 0.055, 0.050 from 
right to left. The results are from DMFT solved in the IPT 
approximation. 

unit-cell volume v. The first three terms arise from ex- 
panding the free energy due to other degrees of freedom 
about a reference cell volume vq , Pq and Bq being the cor- 
responding (reference) pressure and the bulk elastic mod- 
ulus. Expression can be derived from a microscopic 

hamiltonian H = H lat \Ri] - YUja ~ ^j) d \a d o<y + 
U ntfUii, involving both the ion positions Ri and the 
electronic degrees of freedom, when all phonon excita- 
tions are neglected, i.e all lattice displacements R^ — Rj 
are taken to be uniform. It is conventional (see e.g |l6jl 
to use an exponential parametrization for D(v), which 
we linearize since relative changes in v are small: D(v) = 
D exp [-j(v - «o)/uo] — Do [1 - j(v - v )/v ]. As a re- 
sult, the pressure P = —dF/dv and the "compressibility" 
K = -(vdP/dv)" 1 — (vd 2 F/dv 2 )~ 1 are given by: 

P = P - B (v - v )/v Q - { 1 D Q /v )T el (2) 
(K v)- 1 = B /v - (-/Do/vo) 2 Xei (3) 

Here T e i is the (dimcnsionlcss) electronic kinetic energy 
T e i (T,D(v),U) = —dF e i/dD, and Xei is the electronic 
response function: Xei (T, D(v), U) = -d 2 F el /dD 2 . 
Both quantities are associated with the purely electronic 
Hubbard model. Within DMFT, X ei is found 0, @ to 
reach a peak value x™ ax (T, U) for a given temperature 
T at a specific value of D = D^(T, U) , and eventually, 
to diverge at T = Tf, D = Df . This is illustrated in 
Fig.l using DMFT together with iterated perturbation 
theory (IPT) @ in which case Tf ~ 0.02 U, Df ~ 0.4 U. 
In the compressible model, as pointed out in Ref.0, the 
Mott critical endpoint will therefore occur at T c > Tf 
and will be signalled by the divergence of K. From 
this happens when x™i ax nas reached a large enough 
(but not divergent) value such that D x™ ax (T c ) — 
BqVq I (7 2 -Do) . The corresponding critical half-bandwidth 
is D c — Djj"(T c , U) and the critical cell volume is deter- 
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FIG. 2: Temperature dependence of the sound velocity at 
various pressures for parameters corresponding to OC. Inset: 
position and amplitude of anomaly below P c 



mined from (v c — vq)/vq = —{D c — Dq)/('jDq). The 
critical pressure can then be determined using Eq. J2J. 
If T c is close enough to Tf one can use the form (for 
T > Tf): xT x - a el {U)Tf/{T - Tf), which fol- 
lows from mean- field theory |§|. Here, a e i(U) = a/U 
with a ~ 0.5 in our calculations. Hence one obtains: 
AT C /T C = (T C -Tf)/Tf = a[-t 2 D /(B v )](D /U). 

The inverse "compressibility" K~ 1 is directly propor- 
tional to the square of the sound velocity s: s oc 1/y/K. 
At T c , since K diverges, the sound- velocity vanishes. 
Hence, right at the critical point the acoustic phonon 
branch under consideration disperses anomalously. If the 
crystal has inversion symmetry, the dispersion should go 
as lo oc q 2 . Thus from a calculation of the inverse com- 
pressibility along lines similar to that in Ref. one can 
determine the dependence of the sound velocity on tem- 
perature and pressure, including its critical anomaly, as 
reported below. Note that a similar critical behaviour 
of the sound- velocity is well known in the context of the 
usual liquid-gas transition (it| . 

We note that the "pressure" and inverse "compressibil- 
ity" in Eqs. (|2I3I) are associated with the stiffness with 
respect to changes in the unit-cell volume, i.e. to bulk 
strain at fixed number of ions and electrons per unit cell. 
General thermodynamics relates p 2 d/j,/dp (with p the to- 
tal density) to the inverse compressibility provided the 
latter is defined from volume changes at zero strain (i.e 
stemming from vacancy diffusion). Hence, this relation 
does not apply to K (which is related to the stress-stress 
correlation function rather than to the density-density 
one). This also implies that K is not related to the 
"charge compressibility" n e i — dn e i/dp e i (studied the- 
oretically, eg. in (isj). as the latter involves changes 
of the electron density under conditions of zero strain. 
The correspondence between the sound velocity and the 
"compressibility" K as calculated above can be shown to 
be exact for the simplest compressible Hubbard model of 
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FIG. 3: Temperature dependence of the sound velocity 
at various pressures for parameters corresponding to pure 
V2O3. "Negative" pressures can be reached by chromium- 
substitutions: for (Vi-zCr^Os , experiments establish that 
x = 1% corresponds to AP = — 4kbar (~ P c ) (see e.g 0). 



Parameter 


V2O3 


OC 


D 


1 eV 


.13 eV 


vo 


100 A 3 


1700 i 3 


Bo 


2140 kbar 


122 kbar 


7 


3 


5 


B v 


133 eV 


129 eV 


U/D 


2.468 


2.492 


BoV(7 2 A)) 


14.7 


40 



TABLE I: Table of parameter values for 2o3 and the organic- 
conductor k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 (OC) 

phonons modulating the nearest neighbour hopping on 
a cubic lattice, and for longitudinal sound propagation 
along the [111] direction. We hence present it as a rea- 
sonable zeroth- order description of real systems. A more 
realistic description of anisotropic materials should take 
into account the dependence on both polarization and 
direction of propagation. 

The details of the calculational procedure have been 
described in Ref. Q, and we avoid repeating them here. 
The parameter values we have chosen for the two mate- 
rials are given in Table 1. Our approach is to consider 
pure 2o3 and the OC at ambient pressure (i.e., Pq = 1 
bar is essentially zero) as reference compounds, for which 
vq and Bo are taken from experiments. For the OC, we 
take the values Do ~ 0.13 eV and 7 = 5, which were 
found in to be consistent with transport data. For 
2o3, we take Do = 1 eV and 7 = 3 (values considered 
standard for d-electron systems 01). Finally, the value 
of U /Do is adjusted so that the critical pressure is repro- 
duced correctly (P c ~ —4 kbar for 2o3, P c ~ +200 bar 
for the OC). This requires U/Dq to be poised very close 
to the Mott critical value for the pure electronic prob- 
lem (1.26 in our calculations). Using the parameters in 



FIG. 4: The sound velocity as a function of pressure at various 
fixed temperature for parameters corresponding to OC 

Table. 1, one sees that the (dimension-less) combination 
BoVo/{"1 2 D ) is large for both compounds (~ 14.7 for 
2o3 and ~ 40 for the OC). This implies that the Mott 
transition arises for large values of D Xeii i.e. , the ex- 
perimentally observed transition is definitely driven by 
the electronic degrees of freedom, and very close to the 
purely electronic Mott transition. Specifically, we find 
that the relative shift AT C /T C of the critical tempera- 
ture due to the coupling to the lattice, is as small as 
1.4% for 2o3 and even smaller for the organic conduc- 
tors. In fact it is remarkable that despite the very dif- 
ferent values of the bulk modulus of the two materials 
(the OC being much softer than 2o3), the combination 
£>o«o/(7 2 A)) only differs by a factor of 3 between them. 
As a result, the order of magnitude of the sound- velocity 
anomaly is expected to be similar in both materials, since 
K v /Kv = 1 — (j 2 Do/B vo) DoXei- We emphasize that 
there is a significant difference between our choice of pa- 
rameters and that made in Ref. for 2o3: there, a much 
smaller value of Bo (much smaller than the measured ex- 
perimental value) and somewhat larger values of U /D 
were used in order to obtain volume jumps comparable 
to what is seen experimentally, resulting in a large rela- 
tive shift AT C /T C , on the scale of 40%. Because the con- 
tribution of the electronic degrees of freedom to the total 
bulk modulus is comparatively small, it is hard to recon- 
cile the choice made in Ref. ..7] with experimental data 
for this quantity. On the other hand, with our choice of 
Bo, the calculated fractional volume jump Av/vo through 
the transition is ~ 0.2%, too small compared to the ob- 
served jump (~ 1%). We comment on a possible res- 
olution of this problem towards the end of this paper. 

Figs. 2 and 3 show our results for the temperature- 
dependence of the sound velocity computed using the 
method described above, for the parameters represen- 
tative of the OC and 2o3 and for several pressures. At 
the critical pressure P = P c , the sound- velocity vanishes 
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FIG. 5: The sound velocity as a function of pressure at various 
fixed temperatures for parameters corresponding to V2O3. In- 
set: The corresponding phase diagram. The central line and 
the lines on either side correspond to the transition pressure 
Ptr and the spinodal pressures respectively 



according to the mean- field law s oc (T — T^) 1 / 2 (K di- 
verges as 1/(T — T c ), as follows from Eq. J3J). A pro- 
nounced dip remains visible in a rather extended range 
of pressure both above and below P c . The overall shape 
of these curves, as well as the typical order of magnitude 
of the effect are in quite good agreement with the exper- 
imental data for the OC, recently published in For 
example, we obtain a dip of relative size As/sq ~ 10% for 
P ~ 1.3P C , consistent with the experimental observation. 

In Fig. 4 and 5, we show the sound- velocity as a func- 
tion of pressure for various temperatures. This has been 
studied in a less systematic manner in the experiments 
on the OC, but the overall shape and magnitude of our 
results are again consistent with the published data . 
In particular, the curves in Figs. 4-5 show a marked asym- 
metry: on the low-pressure (insulating) side the pressure 
dependence is rather weak and a dip appears only very 
close to P c , while a more gradual and sizeable pressure- 
dependence is observed on the high- pressure (metallic) 
side. This reflects the asymmetry in the electronic re- 
sponse function Xei observed on Fig.l. It is expected 
qualitatively, since electrons participate much less to the 
cohesive energy on the insulating side, and is also ob- 
served experimentally fl3| . 

The inset of Fig. 5 shows the temperature-pressure 
phase diagram for parameters corresponding to 2o3. The 
slope dPtr/dT obtained is too large compared to the ob- 
served value. We believe that this, as well as the prob- 
lem mentioned above, namely the smallness of the calcu- 
lated Av/vq, can be resolved as follows. Let Df r (T) = 
D(vf r {T)) be the half-bandwidth for metal- insulator co- 
existence in the purely electronic problem at tempera- 
ture T. For T well below T c , the fractional differential 
volumes Svi tm /vQ = (vi^ m — vf r (T))/vo of the coexisting 
insulating and metallic phases and the transition pres- 



sure Ptr in the presence of lattice coupling are approxi- 
mately 5vi/v ^ Sv m /v ~ j(T m -7l)/(2B v /D ) and 
P tr ~ B (Df r (T) - £>o)/(7A)) - (lD /v )(T m + %)/2 
(compare |5J| . Here % m are T e i evaluated in the coex- 
isting insulating and metallic phases respectively. These 
expressions are in good agreement with the numerical 
results reported above and in Ref. 0. In the context 
of more realistic models for the electronic problem, one 
expects that differential screening effects would reduce 
(enhance) the effective Hubbard U and correspondingly 
enhance (reduce) T , in the metallic (insulating) phase, 
leading to an overall enhancement of (vi — v m )JvQ. An 
additional enhancement factor might come from the or- 
bital degeneracy in 2o3. These effects could be addressed 
in future studies of multi-band Hubbard models, includ- 
ing screening effects in a self-consistent manner. High- 
precision experimental measurements of optical spectral 
weights in the coexisting insulating and metallic phases 
would also be very interesting, and could provide infor- 
mation on the kinetic energies 7i j771 . 
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